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ABSTRACT 

We present a catalogue of 39 multiple-mergers found using the mergers catalogue 
of the Galaxy Zoo project for z < 0.1 and compare them to corresponding semi- 
analytical galaxies from the Millennium Simulation. We estimate the (volume-limited) 
multi-merger fraction of the local Universe using our sample and find it to be at 
least two orders of magnitude less than binary-mergers - in good agreement with 
the simulations (especially the Munich group). We then investigate the properties 
of galaxies in binary- and multi-mergers (morphologies, colours, stellar masses and 
environment) and compare these results with those predicted by the semi-analytical 
galaxies. We find that multi-mergers favour galaxies with properties typical of elliptical 
morphologies and that this is in qualitative agreement with the models. Studies of 
multi-mergers thus provide an independent (and largely corroborating) test of the 
Millennium semi-analytical models. 

Key words: catalogues - Galaxyfinteractions - galaxies: evolution - galaxies: general 
- galaxies:elliptical and lenticular, cD - galaxies:spiral 



1 INTRODUCTION 



The exact extent to which mergers are able to ac- 
count for the observed properties and morphologies of 
galaxies in the Universe pre occupies much of modern re- 
search fsee lDarg et al . 2010a hereafter DlOa and lDarg et al.l 
l2010bl hereafter DlOb and references therein). If merg- 
ers are not sufficient to explain all observations pertain- 
ing to mass- assembly, then other modes of galaxy for- 
mation must be conjoined to the standard hierarchical 
scheme. If, however, mergers alone (implemented via a 
model of structure formation such as ACDM) are enough 
to explain what we observe in galaxies, then such addi- 
tional processes can (and therefore should) be discarded. 
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In either ca se, mergers are known to occu r at signif- 



icant r ates |Conselice. Raigor fc Mversl 120081: lLotz et al.l 



2008al: iPatton fc Atfieldl |200S|; iBertone fc Conselicd 120091 : 
Stewart et alj|2009al : iRobaina et alJl201Ch and are an indis- 



pensable explanatory resource with respect to galaxy evolu- 
tion. 

To this end the Galaxy Zoo project has established the 
largest homogeneous catalogue of merging galaxies in the lo- 
cal Universe with 3003 visually confirmed mergers (DlOa). 
In this paper we examine an interesting subset of this merg- 
ers catalogue - 39 systems with three or more galaxies merg- 
ing simultaneously (shown in Figure [1} - and use it to es- 
timate the (major) multi-merger fraction of galaxies in the 
local Universe. This fraction and the examined properties 
of these multi-merging galaxies thus provide a useful inde- 
pendent test of semi-analytic models (which are currently 
indispensable aids to the study of galaxy evolution) from 
the typical observables that such models are tuned to repro- 
duce. 

In the ACDM cosmology galactic mergers are 
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Figure 1. Images of the 39 multi-merger systems obtained from GZM1. Each tile has been scaled for optimal viewing, typically with 
sides of ~ 50". 



primarily driven by the coalescence of Dark Mat- 



ter (DM) structures fsee e.e. IWhite & Reesl 


197S; 


Bond et al. 199l|; White & Frenk 1991 


; Lacev & Cold 


1993; 


Jenkins et al.l l200ll; ISDringel et al.l 


2005; IStewart et al. 


2009bh. Once DM haloes have become virialized. fur- 



ther gro wth can only occur through accretion and 
merging dFakhouri fc Mai 120081; iNeistein fc Dekell l20Qg| ; 
iFakhouri. Ma fc Bovlan-Kolchinl l2010t l. The merger histo- 
ries of DM halos can be reproduced through N —Body 
simulations JSpringel et all 120051 ; |Pe Lucia et al.l l2004bl ; 
lHarker et al.l I2006I 1 1 or through analytic appr oximations 



such as P ress-Schechter (IPress fc SchechteJ 19741 ) and its ex 
tensions (|Bond et al.lll99ll ; lBowej|l99ll '). 

To empirically test the accuracy of DM structure for- 
mation, semi-analytic recipes for the evolution of the vis- 
ible matter within these haloes is required. These aim to 
capture the macro-physical processes affecti ng observable 
quan tities such as the Luminosity Function |Benson et al.l 
bOOot ). Each physical process typically involves one or two 
free parameters and so semi-analytic models (SAMs) require 
fine-tuning to reproduce empirical ob servations. The ma - 
jor considerations are p hotoionisation ([Benson et al.ll2002T ). 
shock heating of gas (jCattaneo et al.l 120061) . gas cooling 
dDe Lucia et all |201Ch . AGN feedback (|Bower et all l200rj ; 
ICroton et all l200rj ~ supernovae feedback (also enriching 



1 While the ./V-Body simulations are more or less 'exact' the 
scheme for grouping DM-particles into halos and sub-halos is ar- 
bitrary to a certain extent, in other words, the merger history 
depends on the detailed definition of a halo. 



the IGM with meta l; |Pe Lucia et al.l l2004al ) and mergers 
(|Springel et aljfeOQll ). The resultant 'galaxies' that occupy 
the DM haloes can then be converted to observable lumi- 
nosities through stellar-population synthesis models (e .g. 
iBruzual fc Chariot! bOQSl ; lMarastonlll998l ; lMarastonll2005l 'l. 

Given enough adjustment of the free-parameter val- 
ues and phenomenological ingredients to represent time- 



empirically determined Luminosity Function dCole et al. 


2001; 


Norberg et al.l 200l Huang, et al. 


120031; iPanter et al. 


2004; 


Jones et al.l 120061; iDevereux et al 


l2009f), can be aD- 



proximated arbitrarily well in principle. A major test of 
a model's veracity therefore rests in its capacity to re- 
produce observables that were not involved in it's orig- 
inal calibration. One such test is to determine how 
well a SAM agrees with observed merger rates (or frac- 
tions). Several studies have compared the evolution of 
the galaxy merger rate obtained by the close-pairs tech- 
nique (see i\2.1\ for description) to that of SAMs imple- 
mented in the Millennium Run llKitzbichler fc Whitell2008l ; 
Patton fc Atfieldll2008l; iMateusI I2OO8I ; iHopkins et al.ll201oT ). 



^ — — , , — — 17 ■ — _i_ — - — - — ■ ■ — - — — 1/ 

Bertone fc Conselicell2009 | more recent ly compared the SAM 



of lBertone, De Lucia fc Thomasf 2007 to the merger rate ob 
taine d by the CAS (c o ncentration, asymmet ry and dumpi- 
ness; IConselicd 120031 : IConselice et al] 120031 ) method and 
found the two were roughly consistent for z < 2. 

In this paper, we introduce a new test for 
galaxy-evolution models that has, until now, been 
too difficult to find observationally: the fraction and 
properties of galaxies in (near simultaneous) multi- 
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mergers (mergers of three or more galaxies of sim- 
ilar mass). Severa l individual multi-merger systems 
have been studied dCui et all l200ll ; lAmram et all |2007| ; 
iRines. Finn fc Vikhlininl 120071 ) and a few numerical simu- 
lations of multi-mergers carried out Weil fc Herauistlfl996l : 



lBekkill200ll ; lRenaud. Appleton fc Xu 



20101 ) but no practical 



method has been obtained till now that might locate multi- 
mergers in a near complete manner. 

In Sj2] we describe how the Galaxy Zoo project was 
able to construct such a catalogue in order to estimate 
the multi-merger fraction (carried out in 94.1[l . We then 
compare the binary and multi-merger fractions obtained by 
the Galaxy Zoo project with those of the SAMs of the 
Millennium. They fall broadly into tw o families of mod- 



De Lucia et alJ 120061: 


Bertone. De Lucia & Thomasl 120071) 


and those of Durham 


Cole et al]|2000l:lBenson et al.ll2002l: 


Benson et all 20031; iBaugh et all 120051; iBower et all 120061: 



veloped by both groups and we use those which are pub- 
licly available n For t he MPA, these are the models of 
|Pe Lucia et al]T2006l fc |Pe Lucia fc Blaizotl 120071 (he reafter 



delucia06) and Bertone. De Lucia fc Thomasl 20071 (here 



after bertone 07, an extension of delucia06) and for Durham, 
the model of lBower et al. I l2006l (hereafter bower06). We de- 
scribe some of their characteristics in Ej3l 

By comparing the SAM merger fractions to SDSS ob- 
servations made herein (HJ, we effectively test the accu- 
racy of the build up of dumpiness in the Universe since 
the main factor affecting merger rates is environment. Other 
properties correlate with environment such as stellar mass, 
colour and morphology and so we examine these in both the 
SAMs and Galaxy Zoo catalogues for multi-mergers, binary- 
mergers and single galaxy systems fJJSJ. We summarise our 
results in Sj6] 



2 THE MULTI-MERGER CATALOGUE 
2.1 Finding Multi-Mergers 

As discussed in DlOa, finding mergers amongst surveys 
with ~ 10 6 galaxies is highly non-trivial. Non-parametric 
techniques such as CAS and GM20 (Gini coefficient 
and the second-or der moment of the brightest 20% of 
the galaxy's light; lAbraham. van den Bergh fc Nairl [20031 ; 
lLotz. Primack fc Madaul 120041 ; lLotz et al.ll2008al ) that aim 
to identify parameter spaces uniquely occupied by mergers 
have thus far proved challenging and the prospect of finding 
even more spe cific sub-spa ces limited solely to multi-mergers 
is unrealistic (|Liskerll2008l ; DlOa). 

Likewise, modifying the close-pairs technique (locating 
galaxy pairs within a certain angular separation and redshift 
difference) to find multi-mergers within SDSS is impracti- 
cal due to fibre overlaps. The apparatus within SDSS will 
not acquire spectra (which are needed in such 'blind' meth- 
ods to avoid projection effects) from objects within 55" of 
each other on a single viewing. The conventional close-pairs 



Sec http:/ /www. mpa-garching.mpg.de/millennium/ 




Figure 2. Examples of possible distributions of SDSS objects 
with spectra (red squares) in merging systems. The close-pairs 
technique in particular would be unable to distinguish which of 
these systems is a 'multi-merger' (for example, the top two panels 
have three spectral objects each but the the left-hand system is 
a binary-merger while the right-hand system is a triple-merger; 
the system on the lower panel is a binary-merger but has four 
spectral objects). 



technique is therefore only useful within tile-overlap regions 
jStrauss et al.ll2002l ; lBTanton et alj|2003al ; DlOa). In order to 
find systems with three or more galaxies in the merger stage 
which are within 55" of each other, the system would have 
to rest in part of the sky where there is a double tile-overlap. 

We investigated the practicality of using a modified 
close-pairs technique to find multi-mergers using a close- 
pairs catalogue restricted to 0.005 < 2 < 0.1 that asso- 
ciates spectral galaxy objects with a redshift tolerance of 
0.0017 (co rresponding to a ve locity difference of 500kms~ 
as used in iPatton et ail 120021 ) and projected separation of 
30 kpc. We found that of the 4880 spectral objects in this 
catalogue, only 48 systems (148 objects) had three or more 
spectral objects. Only one of these 48 systems belonged to 
our catalogue indicating how incomplete this technique is for 
finding multi-mergers. Even when three spectral objects do 
fall within these limits, the close-pairs technique would not 
distinguish multi-mergers from normal binary-mergers that 
have been deblended with multiple spectral objects. Figure 
[2] shows examples of this problem - the upper images both 
have three spectral objects, but the left image is a binary- 
merger (with two spectral objects on one galaxy) and the 
right image is a triple-merger. Likewise, without visual in- 
spection, it would be almost impossible to know if the lower 
image, which has four spectral objects, contained two, three 
or four galaxies interacting with each other. 

The only way to find multi-merging systems is therefore 
through visual inspection since humans can readily distin- 
guish between features like bulges and tidal tails and can 
often, on this basis, tell straight away if a system is a multi- 
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Table 1. Catalogue of multi-mergers. 'SDSS Objid' is the unique label for the primary (object 
chosen by Galaxy Zoo) for each system. 'Type' decribes whether the multi-merger is major, middle 
or minor (see text for definition). Mi, M2, etc. give the stellar mass estimates of the galaxies in 
descending order for the system. Galaxies marked * arc part of a major triple- merger (M1/M2 < 3 
and M2/M3 < 3) and are bright enough to be in the volume-limited analysis (M r < —20.55) 
thereby contributing to the major triple- merger fraction (see The entry marked ! has no 
photometric readings and shows a rare failure of the SDSS dcblending routine. 



Type SDSS Objid Ml M2 M3 M4 M5 Ml /M2 M2 /M3 

(x G 1O X M ) 



1 


Minor 


587725491062571042 


11.6 


2 


Major 


587725816951865440 


11.6* 


3 


Major 


587726014546772187 


11.7* 


4 


Minor 


587726016157384716 


11.0 


5 


Major 


587726033303765102 


11.7* 


6 


Middle 


587726033846009923 


11.7 


7 


Minor 


587727865644974418 


11.3 


8 


Minor 


587729385547038748 


12.0 


9 


Minor 


587729772070633570 


11.6 


10 


Minor 


587731511544905794 


11.9 


11 


Major 


587732050018107546 


11.1* 


12 


Major 


587732483809345720 


11.5* 


13 


Minor 


587732483811770549 


11.7 


14 


Minor 


587733081346605098 


10.9 


15 


Major 


587734303805604056 


11.7* 


16 


Minor 


587735696979656747 


10.7 


17 


Minor 


587736477053681696 


11.5 


18 


Major 


587736807771078935 


11.7* 


19 


Major 


587736919972643151 


11.9* 


20 


Middle 


587738952027734025 


11.6 


21 


Major 


587739096450727993 


9.98 


22 


Middle 


587739382058909818 


11.5 


23 


Middle 


587739607547904036 


11.2 


24 


Minor 


587739646203461677 


11.8 


25 


Major 


587739652645191695 


10.4 


26 


Middle 


587739811038101752 


11.4 


27 


Minor 


587739844321542224 


11.4 


28 


Major 


587742014603985009 


11.3* 


29 


Middle 


587742060517064897 


11.9 


30 


Major 


587742062126235884 


11.1* 


31 


Minor 


587744727686381602 


11.3 


32 


Middle 


587745244697329759 


11.5 


33 


Major 


588010359624826929 


11.6* 


34 


Minor 


588013382728221044 


11.5 


35 


Middle 


588013384351678489 


11.1 


36 


Major 


588016878292631762 


11.4* 


37 


Minor 


588017705070100599 


11.4 


38 


Major 


588018090007003216 


11.6* 


39 


Major 


588848898849439845 


11.0* 



merger. The disadvantage with visual identification by indi- 
viduals is, firstly, it is time consuming in general and wholly 
impractical for surveys with ~ 10 6 galaxies and, secondly, it 
is subjective. The Galaxy Zoo project helps overcome both 
of these problems though. It has been shown that human 
classifications averaged over large numbers of individuals 
do provide a good measure of morphol ogical classification 
|Lintott et all |200St iLintott et all |201Ch and merger iden- 
tification (DlOa) and, with instruction, can often identify 
multi-mergers with ease. 

The only method one can use at present is therefore 
to visually flag any system that could be a multi-merger 
(three galaxies of comparable size merging simultaneously), 



10.7 


9.94 






8.06 


5.94 


11.5* 


11.0* 






1.20 


2.65 


11.6* 


11.5* 


11.0 




1.26 


1.51 


10.4 


10.3 






4.21 


1.10 


11.4* 


10.9* 






1.98 


2.88 


11.4 


10.5 






2.17 


7.38 


10.7 


10.1 






3.51 


3.99 


11.1 


10.9 


10.2 




7.56 


1.58 


10.6 


10.5 






8.96 


1.44 


10.9 


! 






8.42 


NA 


10.7 


10.6* 






2.66 


1.08 


11.4* 


11.2* 






1.31 


1.84 


10.7 


10.6 


9.99 




10.4 


1.23 


10.0 


9.18 






8.47 


6.71 


11.4* 


11.3* 






2.01 


1.36 


10.2 


9.62 






3.18 


4.58 


10.7 


10.4 






6.79 


2.00 


11.3* 


11.2* 






2.28 


1.26 


11.7* 


11.2* 


10.6* 




1.93 


2.60 


11.2 


10.1 






2.48 


11.9 


9.83 


9.48 






1.40 


2.24 


11.2 


10.5 






2.14 


4.57 


11.2 


9.90 






1.12 


20.2 


10.8 


10.0 






9.95 


6.71 


10.2 


10.2 


9.38 


8.58 


1.62 


1.02 


11.1 


10.3 






2.33 


6.03 


10.8 


10.8 


10.6 




3.62 


1.03 


11.1* 


11.1* 






1.44 


1.08 


11.4 


10.9 






2.59 


3.16 


10.9 


10.4 






1.80 


2.75 


10.8 


10.2 






3.54 


3.68 


11.4 


10.4 






1.18 


11.0 


11.3* 


11.0* 






2.15 


2.01 


10.9 


10.6 






4.31 


2.18 


11.0 


10.3 


10.0 




1.14 


5.60 


11.0* 


10.6 






2.39 


2.65 


10.7 


10.6 






5.15 


1.19 


11.3* 


11.1 






1.90 


1.58 


10.7* 


10.2* 






1.74 


2.98 



identify the relevant photometry for each galaxy, and use 
this information to more objectively assess what subset of 
those originally flagged are genuinely of comparable size. 
This is explained in more detail in what follows. 
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2.2 Construction of Catalogue 

The measure relevant to mergers produced by the original 
Galaxy Zoo interfac^fl is called the 'weighted-merger-vote 
fraction' (/ m £ [0,1]; DlOa) and was used to create a cat- 
alogue of 3003 mergers called GZM1 (with the constraints 
0.005 < z < 0.1 and 0.4 < f m ). It is calculated by dividing 
the total number of times users classified the system as a 
merger multiplied by a weighting factor representing the re- 
liability of the particular users who examined the system (a 
user is weighted highly if that person tends to agree with the 
majority opinion) and divided by the to tal number of clas- 
sifications that the system received fsee lLintott et all 120081 
for details). 

The creation of this catalogue required the visual re- 
examination of all 3003 systems to check for misclassifica- 
tions and to assign morphologies to the individual galaxies. 
During this process, any system that appeared as though 
it might be a multi-merger was noted for future study. 78 
such systems were flagged and form the parent sample for 
this catalogue. Closer examination of these 78 led to the 
conclusion that several systems were almost certainly not 
merging (but were projection effects), some were too diffi- 
cult to distinguish and 39 of the original 3003 systems, which 
make up the catalogue, are confidently multi-mergers with 
signs of interaction - most of which are discernible in Figure 
m (simple inspection of this figure should convey just how 
variable multi-mergers are in appearance and therefore how 
challenging a pattern-recognition problem it would be to de- 
sign an automated system that could reliably filter these out 
as such). 

Having visually determined that these 39 systems of 
three or more galaxies (with at least one being a Galaxy 
Zoo spectral object associated with an f m value > 0.4) are 
multi-mergers, we manually selected the two (or more) best 
neighbour objects available to represent the other galaxies. 
The 'best' object was judged according to (i) brightness in 
the r-band and (ii) visual common sense (similar to DlOa 
§2.1). Having spectra for at least one galaxy in each multi- 
merging system and photometric data for each individual 
galaxy allows us to calculate rest-frame colours and stellar- 
mass approximations. 



2.3 Mass and Rest-Frame Photometry 
Calculations 

To find a volume-limited major multi-merger fraction we 
need to calculate the rest-frame photometries and stellar 
masses of the galaxies in our catalogue and define what we 
mean by 'major' in this context. 

Stellar masses are calculated as in DlOa by fitting 
the SDSS photometries of each individual galaxy object to 
a library of ste llar-synthesis populations (|Marastonl Il998l; 
iMarastonl [2005ft out to the redshift given by the spectral 
object in the merging system. The photometric errors given 
by SDSS are carried through to estimate errors for the stel- 
lar masses. K-corrected rest-frame colours are calculated as 



3 Note: as of February 2009, a new interface (Zoo Two) has been 
in operation. These results are derived from DlOa whose catalogue 
was constructed using the data obtained in January 2008. 



in DlOb using the photometries and spectral redshifts in- 
put into the publicly a vailable IDL routine kcorrect 4_1_4 
ijBlanton etal] |2003bh . The rest-frame r-band magnitude 
(Mr) is needed to volume-limit the galaxies for our anal- 
ysis as m DlOa and DIObQ 

For each system we refer to the masses of each galaxy 
by descending order where Mi is the most massive galaxy, 
M2 the second most massive, etc. Following DlOa, we call a 
binary- merger a major binary-merger if M2/M1 < 3, else it 
is a minor binary-merger. Similarly, we categorise a multi- 
merger as a major multi-merger if each galaxy of decreasing 
mass is within one third of the mass of the next most massive 
galaxy (M2/M1 < 3, M 3 /M 2 < 3, etc.). Of the 39 systems, 
only 12 are major triple-mergers by this definition. A middle 
multi-merger is where the two most massive galaxies are ma- 
jor (Ma /Mi < 3) but the next two are minor (M 3 /M 2 > 3) 
and a minor multi-merger is where the first two most mas- 
sive galaxies form a minor binary-merger (M2/M1 > 3). A 
summary of the catalogue masses is given in Table [T] Be- 
fore estimating the major multi-merger fraction of the local 
Universe, we now introduce the Millennium SAMs so that 
we can compare the fractions derived from the simulations 
with the empirically observed fraction (0. 



3 THE MILLENNIUM SIMULATION SAMS 
3.1 Dark Matter 

The Millennium Run is an N-body Dark Matter simulation 
of a ACDM Universe. The original version uses ~ 10 10 par- 
ticles in a cubic region of 500ft -1 Mpc sides in comoving 
coordinates and periodic boundary conditions. The simula- 
tion is based on the cosmological parameters fl m = 0.25, 
Ov = 0.75, n b = 0.045, h = 0.73 and as = 0.9 (see 
ISpringel et al.1 [20051 for details). 64 snapshots of the sim- 
ulation were saved (from z = 127 to z = 0) and, from these, 
merger-trees can be constructed. However, the definition of 
haloes and halo substructure is a matter of convention. The 
MPA model defines a Friend-Of-Friend group (defined by 
particles linked to each other by 0.2 times the mean-particle 
separation) and determines substructure using the algorithm 
SUBFIND separating bound structures from unbound parti- 
cles (determined by their velocity relative to their local po- 
tential). 

The Durham models, following lHarker et aill2006l . carry 
out the same analysis with the addition of extra constraints 
on subhalo-definition designed to avoid tenuous 'bridges' 
linking FOF groups but that will disappear in the (near) fu- 
ture rendering them distinct. The subhaloes of a FOF group 
are distinct haloes in their model if (i) the centre of the sub- 
halo is outside twice the half-mass radius of the main halo 
or, (ii) the subhalo has retained more than 75 per cent of 
the mass since the last output time at which it was an in- 
dependent halo. 

The DM build-up is the most important determinant 

4 Specifically, for an object to be capable of obtaining spec- 
tra in SDSS, it must have a petrosian magnitude of r < 17.77 
which means that an object at 2 = 0.1 (the upper-boundary of 
our volume-limited sample) must have a rest-frame magnitude of 
M r < -20.55. 
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Figure 3. Millennium Simulation merger-tree scheme. For all descendant galaxies at snapshot number (SN) =63 with M r < —20.55, 
the progenitors at SN=60 are found and counted to give the number of galaxies in single-, binary- and multi-merger systems. 



with respect to the merging of the baryonic galaxies. When 
a DM (satellite) halo bound to a more massive (main) halo 
becomes sufficiently disrupted and falls below the 20 parti- 
cle limit for a structure to be defined, the countdown begins 
for its central galaxy to merge with that of the main halo. 
The exact timing depends on the SAM details relating to 
dynamical friction that models the decay of the satellite or- 
bit. 



3.2 Baryons 

It is important that the SAM accurately convert baryons to 
photometries if a realistic merger fraction (or rate) is to be 
obtained since any sensible merger fraction must be volume- 
limited (in redshift and magnitude) to ensure completeness. 
SAM recipes ultimately relate DM haloes to galaxy magni- 
tudes in specific bands and if the prescriptions for galaxy 
evolution result in unrealistic stellar populations then the 
volume-limited fractions will be wrong. 

The Durham and MPA SAMs differ on a number of 
details governing the formation and evolution of galaxies. 
Each new model attempts to capture some extra observa- 
tional feature (see lParrv. Eke fc Fr cnk 2009 for a summary 
of the differences and the original papers for details) . For ex- 
ample, delucia06 concentrates on the formation of brightest 
cluster galaxies whereas bertone07 develops the MPA model 
focussing on the treatment of metallicity production and ex- 
change with the IGM. The model reports improvements in 
the suppression of star-formation in small haloes but at the 
expense of galaxy colour-bimodality. 



3.3 Merger Trees 

The SDSS catalogue is limited to the local Universe (z < 
0.1) and, since we assume that the merger rate changes negli- 
gibly over this interval, we are only interested in local merger 
fractions in this paper. The Millennium Database only out- 
puts halos and galaxies at 64 discrete time steps (referred 
to as 'snapshots' labelled by 'snapshot numbers' (SN) for 
to 63) . A schematic is shown in Figure [3] for merger trees 
with the progenitors at SN=60 related to descendents at 
SN= 63. We approximate that a galaxy system would 'look 
like' a merger, on average, if the progenitors of a descendent 
galaxy are identified as individual galaxies at a progenitor 
SN pro9 with redshift corresponding to a look-back time com- 
parable to the time-scale of merger detectability. 

Of course, the vast majority of mergers will be minor- 
mergers whose rates are difficult to constrain observation- 
ally. Since we are interested in comparing the SAM mergers 
to the SDSS catalogue, we consider only progenitors with 
M r < —20.55 and, as far as merger fractions go, we are only 



interested in major mergers (see U). In choosing an 'aver- 
age' or 'typical' time-scale of detectability we face the diffi- 
culty that they depend on the properties of the galaxies in 
the m e rger and the merger-detection techn ique l|Lotz et al.l 



2008b; lLotz et al 



ample, lLotz et al 



2010ai; lLotz et al.ll2010bl ; DlOb). For ex- 



2008b found a median time-scale of 



t n 



r ~ 0.35 ± 0.15Gyr for the detectability of the close- 
pairs technique whereas gas-rich spiral-mergers can remain 
detectable for muc h longer periods w hen using asymmetry 
techniques (> Gvr; lLotz et al.ll2010bh . If a particular popu- 
lation has a high proportion of spiral galaxies compared to 
another, the relative merger fractio ns will therefore appear 
inflated. iBertone fc Conselicd 12009 quote the merger time- 
scale range 0.4 Gyr - 1 Gyr based upon N-Body simula - 
tions and dynamical-friction calculations l|Conselice|[2006l l. 
We choose SN pro9 = 60 for the fiducial snapshot since it 
corresponds to the look-back time ( z ~ 0.064 <-»• 0.6 Gyr) 
closest to the median of the range of these studies (0.35 — 1 
Gyr). We also study how the merger fraction varies with 
SN pr0 g (see Table [2]) and bear this uncertainty in mind in 
interpreting our results. As illustrated in Figure O we allow 
any merger history between SN pro9 and SNd es = 63. 



4 MULTI-MERGER FRACTION OF THE 
LOCAL UNIVERSE 

4.1 SDSS Multi-Merger Fraction 

In DlOa, it was shown how the data from Galaxy Zoo can 
be used to measure the merger fraction of the local Universe 
by effectively integrating over the distribution of f m . Galaxy 
Zoo data is only available for spectral targets and, since a 
meaningful merger-fraction needs to be volume-limited, only 
those galaxies which had acquired spectra were considered 
in the calculation (not all galaxies that are targeted for spec- 
tra obtain them due to fibre clashes in the SDSS apparatus). 
This approach requires special consideration to the fact that 
only about ~ 30% of the SDSS sky has spectroscopic com- 
pleteness in so-called tile overlap regions (a section of the 
sky tha t has been viewed more than once in ta king spec- 
tra; see lStrauss et alJl2002l : iBlanton et alj|2003al : DlOa). In 
single-tile regions, the widths of the spectral fibres on the 
SDSS instrument do not permit spectral targets with in 55" 
of each other to both acquire spectra. One can correct for 
this effect with a multiplicative factor C ~ 1.5 (see Ap- 
pendix A of DlOa) to give a major merger fraction for the 
local Universe of 1 — 3 x C% for galaxies with M r < —20.55. 

We can estimate an upper limit for the multi-merger 
fraction simply by finding the ratio of volume-limited galax- 
ies in major multi-mergers to volume-limited galaxies in 
binary-mergers in GZM1. This catalogue was constructed 
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Table 2. Summary of merger fractions for the Millennium Run and the SDSS fraction estimated 
in this paper. The merger fraction depends on the snapshot number (SNprog) of the progenitors 
which are recorded at redshifts (z) corresponding to look-back times (At) as given. To contribute 
to the numerator of a fraction, a galaxy must be volume-limited (M r < —20.55) and be part of a 
major merger (M2/M1 < 3, M3/M2 < 3, M4/M3 < 3, etc.). *This estimate is based on too small 
a sample to reach any firm conclusions. 
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with f m > 0.4 and it is plausible to assume that multi- 
mergers are amongst the types of merger more likely to 
be classified as 'merging' (in the Galaxy Zoo interface) 
than simple binaries because multi-mergers generally appear 
quite dramatic, prompting the user to go for the merger but- 
ton. The fraction of multi-mergers in systems with f m > 0.4 
is therefore likely to be greater than the fraction of multi- 
mergers in galaxies for all f m . Therefore, by only considering 
fm > 0.4, we can estimate the upper limit of the multi- 
merger fraction in the nearby Universe. 

When we volume limit the 2 x 3003 galaxies in GZM1 
by the constraint M r < —20.55, we are left with 1634 indi- 
vidual galaxies in major mergers (binary or multi). Of the 39 
multi-mergers we have identified, only 16 are major triple- 
mergers and of these systems, 38/48 have M r < —20.55. 
This gives a fraction of 38/1634 and so we can approximate 
the upper limit of the major triple- merger fraction as < 2%0 
This number might be inflated by no more than ~ 50% if 
one takes into account the few systems from the original 78 
(see H2.ip that might have been multi-mergers, but could 
not be resolved sufficiently to be sure. The multi-merger 
to binary-merger ratio appears to be similar, therefore, to 



5 More formally: the ratio of volume-limited (M r < —20.55) 
galaxies in major triple-mergers to volume-limited galaxies in ma- 
jor binary-mergers is < 2%. 



the binary-merger to single-galaxy ratio calculated in DlOa 
(1.5 — 4.5%). We stress that this is a rough estimate for 
the upper boundary, since we have lost information by only 
considering systems for f m > 0.4; but the general result is 
that the probability of finding a galaxy in a merger of N 
galaxies (for iV = 2, 3) is ~ few percent of the probability 
of finding a galaxy in a system of N — 1 galaxies (for these 
volume-limited conditions M r < —20.55 and z < 0.1). 



Extending this empirical query to systems with N ^ 4 
galaxies suffers from small number statistics. Of the 39 
systems, we visually identify only 7 systems as having 
4 or more galaxies merging at once. But none of these 
are majorquadruple-mergers by the appropriate definition: 
M1/M2 < 3, M2/M3 < 3 and M3/M4 < 3. However, the 
systems on lines 3 and 19 of Table [T] only just miss out on 
this definition by having values of M3/M4 slightly above 3 
(with errors making < 3 possible). If both of these systems 
were to be considered as major quadruple-mergers, then the 
ratio of volume-limited galaxies in major quadruple-mergers 
to galaxies in major triple-mergers would be 7±\/7/38 (with 
±\/7 being the poisson counting error). So while the mea- 
sured quadruple-merger fraction is technically zero, it could 
easily have been as high as ~ 25%. Our sample is there- 
fore too small to give accurate merger-fraction estimates for 
systems with 4 or more galaxies of comparable size. 
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4.2 The Millennium Multi-Merger Fraction 

The galaxy databases for the Millennium Run offer much 
larger samples that allow us to calculate merger fractions 
out to quadruple-major mergers. For example, in the model 
of delucia06, at SNd esc = 63 (z = 0), there are 1,673,590 
galaxies with M r < —20.55. The combined number of pro- 
genitor galaxies at SN pros = 60 is 1,943,561. As explained 
above, we classify each descendant galaxy at SN^esc = 63 
as a 'merger remnant' if it has two or more progenitors at a 
given SN pros whose redshift corresponds to a look-back time 
comparable to a merger time-scale (we in fact vary SN pr - 09 
between 58 — 62). 

Of all these 'mergers,' we further classify them as be- 
ing 'major' if all their progenitor masses are constrained by 
our working definition: M2/M1 < 3, M3/M2 < 3, etc. It is 
clear that our Millennium merger-fractions will depend on 
the choice of SN prog since, the lower SN prog is, the greater 
the number of systems there are in mergers (as the merger 
tree branches out with increasing redshift). We therefore 
calculate the merger fractions for a range of SN pro9 values 
and present them in Table [2] referring to fractions derived 
from SN pro9 = 60 (with look-back time of ~ 0.6 Gyr) as the 
fiducial value for each model. 

For the binary-single fraction, we find that all three 
SAMs produce a merger-fraction within the limits of DlOa 
(1.5 — 4.5%). This is in agreement with similar comparisons 
to close-pairs dKitzbichler fc Whitell2008l; Ipatton fc AtfieL ' 



-1 1 



l2008l ; iMateud I2OO8I ) and CAS l|Bertone fc Conselicei I2008 
for the local Universe. However, the SAMs give slightly 
higher percentages for the ratio of galaxies in triple-merger 
systems compared to binary-mergers. We estimated that this 
number should be no more than ~ 2% but the SAMs have 
at least double that fraction for SN pro9 = 60. Only if we 
use the first Millennium time-step SN pros = 62, reducing 
the merger-detectability time-scale to ~ 0.2 Gyr do we get 
agreement between the MPA models and our calculation. 
This might be reconciled by the fact that multi-mergers are 
dominated by elliptical galaxies which have shorter merger- 
time scales (see [JS]). The Durham model predicts roughly 
twice the number of multi-mergers than do the Munich mod- 
els, so the latter are closer to our observations on multi- 
merger fractions. 

For the quadruple-triple ratio the Galaxy Zoo sample is 
too small to offer useful constraints. Comparing the Durham 
to MPA models shows, however, that the prior has a greater 
multi-merger fraction still (over 100% more). 



5. 1 
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Figure 4. Distributions of stellar mass and environment p g for 
galaxies in multi-mergers and galaxies in binary-mergers. The ar- 
rows point at the mean values for the distributions. The peak at 
logp 9 = —3 is artificial (to avoid logO errors we set p g to 10 -3 if 
zero). 
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Figure 5. Members of multimergers (triangles) overlaying the 
catalogue of binary- mergers (grey). 



5 PROPERTIES OF MULTI-MERGING 
GALAXIES 

5.1 SDSS Multi-Merger Properties 

5.1.1 SDSS Morphologies 

One of four morphological categories was assigned to each 
galaxy in the multi-merger systems: S = spiral, E = el- 
liptical, SU = 'unsure' spiral and EU = 'unsure' elliptical. 
Unsure morphologies are common, especially in late stage, 
distant systems where structural indicators like spiral arms 
cannot be seen. The same four categories were used in DlOa 
where it was found that spirals (S and SU) outnumber ellip- 



ticals (E and EU) in volume limited galaxies in mergers by at 
least 3:1. This is about twice the ratio of the global popula- 
tion and it was argued in DlOb that this is most likely due to 
the fact that mergers involving spirals remain detectable for 
longer periods than mergers involving ellipticals. In multi- 
mergers, by contrast, we found that of the spiral to elliptical 
ratio for volume-limited galaxies was closer to 1:1. In other 
words, multi-mergers appear to favour ellipticals in our sam- 
ple compared with their occurrence in binary-mergers and 
the global galaxy population. 
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5.1.2 SDSS Environment, Colours and Stellar Masses 

It seems unlikely that this high occurrence of ellipticals 
in multi-mergers is entirely due to a selection effect. It is 
possible that bulge dominated galaxies can be identified in 
multi-mergers longer than spirals since the bulge is a key 
feature indicating how many galaxies were originally in- 
volved in a multi-merger. Spirals remain detectable in merg- 
ers for longer (DlOb) so the duration of detectability for 
multi-mergers involving spirals should be longer still than 
for multi-mergers involving ellipticals. Despite this, the fact 
that ellipticals feature so prominently in multi-mergers sug- 
gests that environment has some influence, i.e. dense envi- 
ronments are more likely to host multi-mergers. It was found 
in DlOb that (binary) mergers tend to occupy slightly denser 
environments than galaxies in the global population and so 
it seems that merger number scales with environment. 

We measure the environment of the multi-mergers di- 
rectly using the adaptive Gaussian environment parame- 
ter, p g , following DlOb. This is a sophisticated measure of 
galaxy number density that utilises all the relevant spectral 
information in the SDSS database jSchawinski et all 120071 ; 
DlOb) returning a value p g £ R + for any coordinate in the 
SDSS sky (ra,dec,z). The parameter p g (ra, dec, z, a) starts 
by finding close neighbours within an initi al radius of a for 
each galaxy (we use a — 2.0Mpc following ISchawinski et all 
12007( 1. It then adapts this radius depending on the initial re- 
turn in order to compensate for the "finger- of-God" effect 
and is weighted such that p g increases the nearer its neigh- 
bours are. If p g = then there are no neighbours within a a 
radius. A galaxy with p g = 1 roughly corresponds to it being 
at the centre of a sphere of r adius 3 Mpc with ten g alaxies 
randomly distributed within l)Schawinski et aDl2007f j. 

Figure [4] shows the distribution of p g for the multi- 
merger and binary- merger systems (top- left panel). The 
multi-mergers on average occupy environments with slightly 
higher values of p g . The Kolmogorov-Smirnov statistic be- 
tween the binary and multi-merger data sets provides a mea- 
sure of the difference between their cumulative distributions 
and we find them to be different with a significance level of 
> 99%. 

The observed multi-mergers are therefore found in more 
dense environments though the primary reason is not clear: 
do multi-mergers favour ellipticals because they are likely 
to take place in dense environments or are multi-mergers 
with ellipticals easier to visually identify and therefore make 
multi-mergers appear to favour dense environments? Proba- 
bly both factors are at play and it is difficult at this stage 
to disentangle the effects quantitatively. 

Likewise, Figure 0] indicates that the stellar masses of 
the galaxies in multi-mergers are on average greater than 
their binary-merger counterparts. We compare the SDSS 
masses directly to the Millennium SAMs in £15.2.21 Colour is 
yet another quantity that correlates with morphology and, 
as Figure[S]shows, the u—r colours of multi-merging galaxies 
are redder than their binary-merger counterparts. The em- 
pirical evidence is therefore emphatically clear that galaxies 
observed in multi-mergers are more likely to be early-type 
than single and binary-merger galaxies. 
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Figure 6. Distributions of Bulge- Total mass ratio of galaxies in 
single-, binary- and multi- merger systems for all three SAMs at 
SN prog = 60. The arrows indicate the mean values for each sam- 
ple. All three predict that binary and multi-mergers have more 
'elliptical' like morphologies compared with isolated systems - a 
likely concomitant of the favourability of mergers to occur in high- 
density environments. 



5.2 Millennium Multi-Merger Properties 

5.2.1 Millennium Morphologies 

The morphologies of the SDSS multi-mergers and the Mil- 
lennium multi-mergers cannot be compared directly since 
the SDSS morphologies are obtained visually. However, 
we can determine the qualitative relationship of the Mil- 
lennium galaxy morphologies using the Bulge- Total stel- 
lar mass ratio as a proxy (with a high ratio correspond- 
ing to ellipticals and low ratio to spi rals). Several studies 
have d efined morph o logies t his way jKhochfar fc Burkert 
20031: iBenson etal] 120071; iBertone. De Lucia fc Thomas 



20071 : iParrv. Eke fc Frenkl 120091 ) ■ IBenson et alj 120071 found 
that the SDSS and Durham SAM produced qualitatively 
similar disc-bulge Luminosity Functions (see their Figure 
17). The distributions of Bulge/Total mass for the galax- 
ies of the three SAMs at SN pr . 09 = 60 is shown in Figure 
[S] All three models produce the expected qualitative result 
that the more galaxies there are in a (merger) system, the 
more bulge dominated they will be. The Durham model has 
systems generally more b ulge-dominated in agreement with 
IParrv. Eke fc Frenldl2009l . 
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As argued in i )5,1.2l binary- and multi-mergers oc- 
cur most favourably in high-density environments and this 
is where interaction s (inducing gravitatio nal torques, see 
iHopkins et"ai1l2009al ; iHopkins etai1l2009bl ) causing disk in- 
stability are common. This is no doubt largely responsible 
for the well established morp hology-environ ment relation- 
ship (established since at least lDressler|[l98(ih and so merg- 
ers, taking place more favourably in denser environments, 
are more likely to be elliptical. 



5.2.2 Millennium Environment, Colour and Stellar 
Masses 

The Millennium SAM catalogues do not provide a direct 
measure of environment and so we decline to test this prop- 
erty, though it seems almost certain that multi-mergers will 
favour high-density environments given the bulge-dominated 
morphologies exhibited by the SAMs as discussed in £15.2,11 

The colours of delucia06 and bert one07 were examined 
in lBertone. De Lucia fc Thomasll2007l and it was found that 
bertone07 did not reproduce the colour bi-modality of the 
local Universe as well as delucia06. The model of bowerOB 
does reproduce th e local bi-modality well (see e.g. Figure 4 
iBower et al.ll200rt ). However, since the colours are derived 
secondarily from the stellar populations comprising their 
putative galaxies (all us ing the stellar synthesis models of 
iBruzual fc CharlotllioO^ ). we analyse their stellar mass dis- 
tributions as a main form of comparison. 

Figure [7] shows t he mass distributions for the SPSS a nd 
Millennium galaxies. iBertone. De Lucia fc Thomasll2o"o"7l re- 
ported that both their model and that of delucia06 slightly 
underestimate the Luminosity Function as a function of stel- 
lar mass in the range 10.6 < logM* < 11.2. This helps ex- 
plain the very slight discrepancy between the single-galaxy 
mass distributions for the MPA models (note the black ver- 
tical arrows show the mean masses of the SAMs are ~ 1 
dex less than that of the SDSS single systems). The top 
and bottom panels show the SDSS distributions where the 
data points on the bottom are allowed to vary along the 
poisson-count error bars in order to most closely match the 
Millennium distributions. Also, importantly, the distribu- 
tions for the binary-mergers, using the catalogue of DlOa, 
are adjusted so as to include a spiral-elliptical ratio of 3 : 2 
since it was argued in DlOb that spirals are over-observed 
by f actor ~ 2 due to l onger time-scales of detectability (see 
also lLotz et alj|2008bl ) . 

By contrast, it was suggested in £15.1.11 that ellipticals 
might be over-observed but since we cannot quantify this, we 
cannot correct for it. It suggests that the multi-merger mass 
distributions might not increase at quite the rate suggested 
by the SDSS distributions. All this considered, the MPA 
models both do well in reproducing the fact (with slight 
underestimation) that galaxies in binary and multi-mergers 
increase in mass by ~ 1 dex per extra galaxy in the system 
(on average). By contrast, the model of bower06 appears 
to reproduce the single-galaxy mass distribution to great 
accuracy but shows very little in the way of increasing mass 
with merging. 



6 SUMMARY 

Through the Galaxy Zoo project, we have assembled a cat- 
alogue of multi-merger systems in the local Universe (z < 
0.1). Multi- mergers can only be found in SDSS through vi- 
sual inspection (currently) since close-pairs suffer from fiber 
overlaps and automated methods like CAS and GM20 are 
not sensitive enough to consign a given image to a 'multi- 
merger' parameter space (distinct from binary-mergers). 
The original Galaxy Zoo interface was not setup specifically 
for the task of finding multi-mergers, but will be in the fu- 
ture. This should produce a larger, more-complete catalogue. 

Nonetheless, we argued that our catalogue is sufficient 
to provide a rough estimate of the (major) multi-merger 
fraction and gave an upper bound such that < 2% of all 
volume-limited galaxies (M r < —20.55) in a major merger 
(most of them being binary) are specifically in a multi- 
merger. This is about the same percentage as that calcu- 
lated for the (binary) major merger fraction (1.5 — 4.5%) 
for the same volume-limiting constraints (DlOa). However, 
our sample is not large enough to find the quadruple-major 
merger fraction (or beyond). 

The Millennium SAMs gave similar merger fractions 
for the number of volume-limited galaxies in major binary- 
mergers compared to single galaxies (~ 2%). However, the 
next level of merger exhibited some disagreement between 
the SAMs and observation with the SAMs over-predicting 
galaxies in multi-mergers by at least factor ~ 2. The Durham 
model offered a multi-merger fraction twice that of the MPA 
model. However, since we have shown that galaxies in multi- 
mergers tend to be elliptical, which have shorter time-scales 
of detectability, this could justify taking the Millennium 
multi-merger step at SN pro9 = 62 in which case the frac- 
tion would be within our rough observational limit in the 
bertone07 model. 

Examining the properties of the galaxies in these multi- 
mergers to those in binary-mergers and single galaxy sys- 
tems, we found that the (volume-limited) galaxies in multi- 
mergers have greater stellar masses (Figure [7J, redder 
colours (Figure [5} and occupied denser environments (Fig- 
ure 2| . Such properties are characteristic of elliptical galax- 
ies and we found a high proportion of ellipticals-spirals (at 
about 1 : 1) in multi-mergers compared to the single sys- 
tems (at about 2 : 3) and binary-mergers (at about 1:3). 
We argued that this is unlikely to be entirely due to a se- 
lection effect. To corroborate this, we compared our results 
with the SAMs. 

We found good qualitative agreement: all three mod- 
els predicted that the bulge-total mass ratio increases with 
merger-number corroborating the observation that multi- 
mergers favour ellipticals and, implicitly, occupy denser envi- 
ronments on average compared to binary-mergers and single 
galaxies. The MPA models also agreed qualitatively with the 
fact that galaxy mass increases, on average, with the num- 
ber of galaxies in the system. The Durham model appeared 
to slightly underestimate this effect. 

Future work with the Galaxy Zoo project aims to ex- 
pand and improve these multi-merger catalogues as they 
provide a valuable, independent test to semi-analytical mod- 
els. 
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Figure 7. Mass distributions for the single-, binary- and multi-merger systems in SDSS and Millennium. The left, middle and right 

columns correspond to the delucia06, bertone07 and bowcr06 models respectively. The top panels show the mass distributions for the 

SDSS systems. The bottom panels show the same (SDSS) data except the plotted points are allowed to vary along the error bars 

(representing the poisson counting errors) in order to minimise a \ 2 significance test with the Millennium model (in the middle panel). 

The arrows show the mean values for each distribution. All galaxies arc volume-limited (z < 0.1, M T < —20.55) and the progenitors are 

taken at SN prog = 60. The Millennium SAM masses are determined here using h~ 1 = 1/0.73. 
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